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(Dated: March 29, 2001) 
We have observed the decays B — > <f>K and (f>K* in a sample of over 45 million B mesons 
collected with the BABAR detector at the PEP-II collider. The measured branching fractions 
are B(B+ -> (j>K+) = (7.7±^ ± 0.8) x 1CT 6 , B{B° -> 0K ) = (S.ltf.l ± 0.8) x lO" 6 , B{B+ -> 
0/C+) = (9.7ia;4 ± 1.7) x 10" 6 , and B(B° -> <?!>A"* ) = (8.6iai ± 1.1) X 10" 6 . We also report 
the upper limit B(B + -> 0tt+) < 1.4 x 10~ 6 (90% CL). 

PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 



The decays of B mesons into charmless hadronic fi- 
nal states provide important information for the study 
of CP violation and the search for new physics. Decays 
into final states containing a <fi meson are particularly 
interesting because they are dominated by b — > s(d)ss 
penguins (Fig. |l|), with gluonic and electroweak contri- 
butions, while other Standard Model contributions are 
highly suppressed jj| . These modes thus provide a direct 
measurement of the penguin process, with potential ben- 
efits to estimates of direct CP violation. They also allow 
an independent measurement of sin 2/3 0). Comparison 
of the value of sin 2/3 obtained from these modes with 
that from charmonium modes, as well as various tests of 
isospin relationships, can probe for new physics |J. 

In this paper we present measurements of four such 
decays: B+ -> <j>K+, B° -► (f)K°, B+ -► (f>K*+ , and 
B° — ► <fiK*°. Charge conjugate states are assumed 
throughout this paper and measured branching fractions 
are averaged accordingly. The <f>K + and <fiK*° modes 
have been previously seen [||. 

The data were collected with the BABAR detector || 
at the PEP-II asymmetric e + e~ collider H located at 
the Stanford Linear Accelerator Center. The results pre- 
sented in this paper are based on data taken in the 1999- 
2000 run. An integrated luminosity of 20.7 fb _1 was 
recorded corresponding to 22.7 million BB pairs at the 
T(4S) resonance ( "on-resonance" ) and 2.6 fb _1 about 
40 MeV below this energy ("off-resonance"). 

The asymmetric beam configuration in the laboratory 
frame provides a boost to the T(4S) increasing the mo- 
mentum range of the i?-meson decay products up to 




u, d u, d it, d u, d 

(a) (b) 



FIG. 1: Gluonic penguin diagrams describing the decays 
B — > <f>K, <j>K* , and <f)n: (a) internal (b) flavor-singlet. 



4.3 GeV/c. Charged particles are detected and their mo- 
menta are measured by a combination of a silicon vertex 
tracker (SVT) consisting of five double-sided layers and 
a 40-layer central drift chamber (DCH), both operating 
in a 1.5 T solenoidal magnetic field. With the SVT, a 
position resolution of about 40 /jm is achieved for the 
highest momentum charged particles near the interac- 
tion point, allowing the precise determination of decay 
vertices. The tracking system covers 92% of the solid an- 
gle in the center-of-mass system (CM). The track finding 
efficiency is, on average, (98±1)% for momenta above 
0.2 GeV/c and polar angle greater than 0.5 rad. Pho- 
tons are detected by a Csl electromagnetic calorimeter 
(EMC), which provides excellent angular and energy res- 
olution with high efficiency for energies above 20 MeV [§. 

Charged particle identification is provided by the aver- 
age energy loss (dE/dx) in the tracking devices and by a 
unique, internally reflecting ring imaging Cherenkov de- 
tector (DIRC) covering the central region. A Cherenkov 
angle K-tt separation of better than 4er is achieved for 
tracks below 3 GeV/c momentum, decreasing to 2.5<r at 
the highest momenta in our final states. Electrons are 
identified with the use of the EMC. 

Hadronic events are selected based on track multiplic- 
ity and event topology. We fully reconstruct B meson 
candidates from their charged and neutral decay prod- 
ucts, where we recover the intermediate states tt° — ► 77, 
K° -> K° s -> 7T+7T-, -> K+K-, K*+ -> K tt+ or 
K + tt°, and K*° — * K + ir~ . Candidate charged tracks are 
required to originate from the interaction point (within 
10 cm along the beam direction and 1.5 cm in the trans- 
verse plane), and to have at least 12 DCH hits and a 
minimum transverse momentum of 0.1 GeV/c. Looser 
criteria are applied to tracks forming Kg candidates to 
allow for displaced decay vertices. Kaon tracks are dis- 
tinguished from pion and proton tracks via a likelihood 
ratio that includes, for momenta below 0.7 GeV/c, dE/dx 
information from the SVT and DCH, and, for higher mo- 
menta, the Cherenkov angle and number of photons as 
measured by the DIRC. A kaon (pion) candidate is any 
track not identified as a proton or pion (kaon) . 

We reconstruct ir° mesons as pairs of photons with 
a minimum energy deposition of 30 MeV. The typical 
width of the reconstructed 7r° mass is 7 MeV/c 2 . A 
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FIG. 2: The two-kaon invariant mass in the <j> signal region 
(left). Superimposed to the data is the fit to a relativistic 
P-wave Breit-Wigner for the <f> convoluted with a Gaussian 
on top of a polynomial background. The mass resolution is 
f.l MeV/c 2 . The plot to the right shows a Breit-Wigner fit 
to the K°tv + invariant mass in the K* + signal region. Both 
fits use Breit-Wigner parameters from Ref. ||. 

±15 MeV/c 2 interval is applied to select tt° candidates. 

We combine pairs of tracks with opposite charge from 
a common vertex to form Kg, tf>, and K*° candidates. 
The selection of K s candidates is based on the invariant 
two-pion mass i\M vv — m K o\ < 10 MeV/c 2 ), the angle 
a between the reconstructed flight and momentum di- 
rections in the plane transverse to the beam direction 
(cos a > 0.999), and the measured lifetime significance 
(t/<t t > 3). For the softer K s from K* + decays we relax 
the criteria to 12 MeV/c 2 and cos a > 0.995. 

For <j) candidates, both daughters are required to be 
kaon candidates. The invariant mass for the K + K~ pair 
must lie within 30 MeV/c 2 of the <ft mass (see Fig. ||). 

The natural width of the K* dominates the resolution 
in the invariant mass spectrum. For K*° candidates the 
Kn invariant mass interval is ±100 MeV/c 2 . The selec- 
tion of K* + comprises K + ir° and Kgir + combinations 
within a Kn mass interval of ±150 MeV/c 2 (see Fig. |). 
We require particle identification for the charged daugh- 
ters of the K*. To suppress combinatorial background wc 
restrict the K* + — > K + tt° helicity angle (cos9 H > — 0.5 
as defined below). This effectively requires the 7r° mo- 
mentum to be above 0.35 GeV/c. 

The helicity angle 9h of a (f> or K* is defined as the 
angle between one resonance daughter direction and the 
parent B direction in the resonance rest frame. For 
pseudoscalar-vector B decay modes, angular momentum 
conservation results in a cos 2 9h distribution, whereas in 
decays into vector- vector states, the distribution is the re- 
sult of an a priori unknown superposition of transverse 
and longitudinal polarizations. 

We identify B meson candidates kinematically us- 
ing two independent variables ||, AE = (EiEs — 
Pi • Pb — s /^)/v^ ano - the energy-substituted mass 
77iEs = y/(s/2 + Pi ■ p B ) 2 /Ef - pf, where y/s is the to- 
tal e + e~ CM energy. The initial state four- momentum 



(Ei,pi) derived from the beam kinematics and the four- 
momentum (Eb,Pb) of the reconstructed B candidate 
are all defined in the laboratory. The calculation of m^g 
only involves the three-momenta of the decay products, 
and is therefore independent of the masses assigned to 
them. For signal events AE peaks at zero and Toes at 
the B mass. Our initial selection requires \AE\ < 0.23 
GeV and m ES > 5.2 GeV/c 2 . 

Monte Carlo (MC) simulation [jl0| demonstrates that 
contamination from other B decays is negligible. How- 
ever, charmless hadronic modes suffer from large back- 
grounds due to random combinations of tracks produced 
in the quark-antiquark (qq) continuum. The distinguish- 
ing feature of such backgrounds is their characteristic 
event shape resulting from the two-jet production mech- 
anism. We have considered a variety of event shape vari- 
ables in the CM that exploit this difference. 

One such variable is the angle 9t between the thrust 
axis of the B candidate and the thrust axis of the rest of 
the event, where the thrust axis is defined as the axis that 
maximizes the sum of the magnitudes of the longitudi- 
nal momenta. This angle is small for continuum events, 
where the J3-candidate daughters tend to lie in the qq 
jets, and uniformly distributed for true BB events. Thus 
we require | cos6» T | < 0.9 (0.8 for 4>K*+). 

Other quantities that characterize the event shape are 
the B polar angle 9b and the angle 9 q q of the £>-candidate 
thrust axis, both defined with respect to the beam axis, 
as well as the angular energy flow of the charged parti- 
cles and photons relative to the £?-candidate thrust axis. 
For T(4S) decays into two pseudoscalar B mesons, the 
9b distribution has a sin 2 9b dependence, whereas the 
jets from continuum events lead to a uniform distribu- 
tion in cos9b- In 9 q q, the continuum jets give rise to a 
(1 + cos 2 9 q q) distribution, while the thrust direction of 
true B decays is random. We enhance the background 
suppression by forming an optimized linear combina- 
tion of eleven variables (Fisher discriminant): |cos#b|, 
|cos#gg|, and energy flow into the nine 10° polar angle 
intervals coaxial around the B candidate thrust axis B . 

We use an extended unbinned maximum likelihood 
(ML) fit to extract signal yields. The extended likeli- 
hood for a sample of N events is 

(M \ N / M \ 
~^2 Ui ) IT ( ^niViixjid) J , (1) 
i=i / j=i \i=i / 

where Vi(xj;a) describes the probability for candidate 
event j to belong to category i, based on its measured 
variables Xj, and fixed parameters a that describe the 
expected distributions of these variables in each of the 
M categories. In the simplest case, the probabilities are 
summed over two categories (M = 2), signal and back- 
ground. The decays B + — > (j>K + and B + — > (j>ir + are 
fit simultaneously with two signal and two correspond- 
ing background categories [M = 4). The event yields rii 
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in each category are obtained by maximizing C [0. Sta- 
tistical errors correspond to unit changes in the quantity 
X 2 = —2 In C around its minimum value. The significance 
of a signal is defined by the square root of the change in 
X 2 when constraining the number of signal events to zero 
in the likelihood fit. 

The probability Vi (xj ; a) for a given event j is the 
product of independent probability density functions 
(PDFs) in each of the fit input variables Xj. These are 
AE, ?TiEs, Mkk for all channels, Mk-k for the <f>K* chan- 
nels, the 4> helicity angle for pseudoscalar-vector decays, 
and event shape quantities as discussed below. For the 
simultaneous fit to the decays B + — > <j>K + and <j)ir + we 
include normalized residuals derived from the difference 
between measured and expected DIRC Cherenkov angles 
for the charged primary daughter. Additional separation 
between the two final states is provided by AE. 

The fixed parameters a describing the PDFs are ex- 
tracted from signal and background distributions from 
MC simulation, on-resonance AE—rriES sidebands, and 
off-resonance data. The MC resolutions are adjusted by 
comparisons of data and simulation in abundant calibra- 
tion channels with similar kinematics and topology, such 
as B — > Dw, Dp with D — > Kir, Kirir. The simulation 
reproduces the event-shape variable distributions found 
in data. The Cherenkov angle residual parameterizations 
are determined from samples of D° — > K~ir + originating 
from D* decays. 

For the parameterization of the PDFs for AE, ttieSj 
and resonance masses we employ Gaussian and Brcit- 
Wigner functions to describe the signal distributions. For 
the background we use low-degree polynomials or, in the 
case of rriES; an empirical phase-space function The 
background parameterizations for Mkk and Mk-k also 
include a resonant component to account for <j) and K* 
production in the continuum. The <pK (4>ir) helicity- 
angle distribution is assumed to be cos 2 Oh for signal. 
The background shape is again separated into contri- 
butions from combinatorics and from real 4> mesons, 
both fit by nearly constant low-degree polynomials. The 
Cherenkov angle residual PDFs are Gaussian for both the 
pion and kaon distributions. The thrust and production 
angle PDFs are parameterized by polynomials, with the 
exception of the background in | cos Ot \ , where we use an 
exponential. The Fisher discriminant is described by an 
asymmetric Gaussian for both signal and background. 

For all modes, we test the fit response for various 
choices of preselection and fit strategies with samples 
generated according to the PDFs, each containing the ex- 
pected number of events in signal and background. Signal 
yields were found to be unbiased. In the (f>K° analysis 
the results of our tests show that fitting either to | cos Ot \ 
and cos Ob or to the Fisher discriminant yields compara- 
ble significance. Thus, we use only the thrust and B 
polar angle in this analysis. In the other modes we find 
that the additional background discrimination provided 



TABLE I: Summary of results; e denotes the reconstruc- 
tion efficiency and £tot the total efficiency including daugh- 
ter branching fractions, both in percent; N is the number of 
events entering the ML fit, n s i g the fitted number of signal 
events, S the statistical significance (in Gaussian a), and B 
the measured branching fraction including statistical and sys- 
tematic errors. The subscripts in the (f>K* + modes refer to 
the kaon daughter of the K* + . 



Mode 


e 


£tot 


N 




S 


6(10~ 6 ) 


(j>K+ 
<f>K° 
<j)K* + 


36.4 
37.4 


17.9 
6.1 
4.9 


4202 
351 


31.4+ 6 5 .9 
10.8±« 


10.5 
6.4 
4.5 


7.71^4 ±0.8 
8.1^5 ±0.8 
9.7tt;t ± 1.7 


<t> K K + + 
<t> K KO 


15.1 
21.5 


2.5 
2.4 


781 
381 


7.1- 4 3' 3 4 

4-4±»:S 


2.7 
3.6 


12.8l^±3.2 
8.0^;? ± 1.3 


<f>K*° 


26.3 
38.9 


8.6 
19.1 


2517 
4202 


16.9+1;? 
9+ 2 1 


6.6 
0.6 


8.6±H ± 1-1 
< 1.4 (90% CL) 



by the Fisher discriminant improves the expected signif- 
icances of the results, and use this approach. 

The results of our ML fit analyses are summarized in 
Table |. For the branching fractions we assume equal 
production rates of B°B° and B + B~ . We find significant 
signals in all four B — > <j>K and 4>K* decay modes. The 
number of fit events, their statistical significance, and 
the ML fit x 2 values are well reproduced with generated 
samples. Projections of the input variables are in good 
agreement with the fit results, as shown for tties m Fig. |^. 

We check the stability of our results by reducing the 
number of input variables in the fit. In particular, we 
find statistically significant signals even if event shape 
variables are omitted from the fit and only preselection 
criteria are required. Correlations among the input vari- 
ables are found to be less than 10%. 

Systematic uncertainties in the ML fit originate from 
assumptions about the signal and background distribu- 
tions. We vary the PDF parameters within their respec- 
tive uncertainties, and derive the associated systematic 
errors. They range between 4 and 9% (17% for the final 
state that includes a 7r°). The signals remain statistically 
significant within these variations. 

The dominant systematic errors in the efficiency are 
track finding (1.2% per track), particle identification 
(2% per track), and Kg and ir° reconstruction (7% and 
5%, respectively). Other minor systematic effects from 
event selection criteria, daughter branching fractions, 
MC statistics, and B meson counting sum to less than 
4%. The efficiency in the ML fit to signal samples can be 
less than 100% because of fake combinations passing the 
selection criteria, and we account for this with a system- 
atic uncertainty (2-5%). This effect is larger in the K* 
final states because of broader distributions and combina- 
torial ir° background. Uncertainties in the efficiency only 
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FIG. 3: Projections onto the variable m^s- The histograms 
show data for (a) B+ -» (b) B° -> 4>K°; (c) B+ -> 

(f>K* + ; (d) B — > <f>K*° after a requirement on the signal 
probability Vmg./SVi with the PDF for mES excluded. In (c) 
the histogram is the sum of the two <f>K* + channels while the 
shaded area is K*+ -> A-°7r + alone. The solid (dashed) line 
shows the PDF projection of the full fit (background only). 



affect the branching fraction, but not the significance of 
a result. 

In the vector-vector final states we average efficien- 
cies for the transverse and longitudinal angular polariza- 
tions and assign a systematic error as the rms spread of a 
uniform efficiency distribution between the two extreme 
cases (6% in 4>K*+ , 14% in 4>K*+ + , and 2% in (f>K*°). 
We combine the results from the two K* + decay channels 
using \ 2 distributions convoluted with the uncorrelated 
part of the systematic errors. 

The fit result for the B + — > <frir + branching fraction 
is (2.1±2 'i ± 0.5) x 10~ 7 . Given the signal yield of less 
than one event, we quote an upper limit obtained by inte- 
grating the normalized likelihood distribution. The limit 
incorporates changes by one standard deviation from un- 
certainties in PDFs and the reconstruction efficiency. 

Event counting analyses, based on the same variable 
set Xj as used in the fits, serve as cross-checks for the ML 
fit results. The variable ranges are generally chosen to 
be tighter in order to optimize the signal-to-background 
ratio, or upper limit, for the expected branching frac- 
tions. We count events in a rectangular signal region in 
the AE-ities plane, and estimate the background from a 
sideband area. For B + — > <fiK + we find 43 events in the 
signal region (expected background 9.4); the correspond- 
ing numbers are 10 (2.8) for (pK°, 6 (2.2) for cf>K*+, 22 
(7.3) for 4>K*°, and 2 (3) for (f>ir+ . The branching frac- 
tions measured using this technique are in good agree- 
ment with those arising from the ML fit analysis. 

In summary, we have observed B decays to <pK + , (j)K°, 
4>K* + , and 4>K*° with significances, including systematic 



uncertainties, of greater than four standard deviations 
(Table ||). The agreement between the branching frac- 
tions of charged and neutral modes is in accordance with 
isospin invariance under the assumption of penguin dia- 
gram dominance. The decay B + — > </>7r + has both CKM 
and color suppression relative to cj)K + B and is therefore 
not expected to be observed in the present data sample. 
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